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Abstract 

We derive a modified form of the BFKL equation which enables the structure of 
the gluon emissions to be studied in small x deep inelastic scattering. The equation 
incorporates the resummation of the virtual and unresolved real gluon emissions. We 
solve the equation to calculate the number of small x deep-inelastic events containing 0, 
1,2... resolved gluon jets, that is jets with transverse momenta qx > \x. We study the 
jet decomposition for different choices of the jet resolution parameter jx. 



1. Introduction 

The advent of the HERA electron-proton collider has opened up the possibility of testing 
QCD in the new and hitherto unexplored small x regime. The HERA measurements of the 
proton structure function F 2 (x, Q 2 ) show a striking rise with decreasing x, which with the latest 
data is now known with considerable precision [[[], [2|. On the other hand from the theoretical 
point of view we know for sufficiently small x, such that a^lnl/x ~ 1, that it is necessary to 
resum the (c^lnl/a;)™ contributions in order to obtain reliable perturbative QCD predictions. 
At leading order this is accomplished by the Balitzkij, Fadin, Kuraev, Lipatov (BFKL) equation 
0. This equation effectively corresponds to the sum of gluon ladder diagrams of the type shown 
in Fig. 1 in which the transverse momenta qr are unordered along the chain. This should be 
contrasted with DGLAP evolution where, in the leading InQ 2 approximation, the transverse 
momenta are strongly ordered from the hadronic to the hard scale Q 2 which in deep-inelastic 
lepton scattering is provided by the virtuality of the photon, namely 

Q 2 > k 2 T > k 2 nT > ... (1) 

Both BFKL and DGLAP evolution lead to an increase of the deep-inelastic scattering structure 
functions with decreasing x. In fact it is possible to obtain a satisfactory description of the rise 
of the structure function measured in the HERA small x regime using both approaches [S]. |(| . 
The inclusive nature of the structure function F 2 makes it extremely difficult, even with the 
precise HERA data, to use the observed x behaviour to reveal the underlying dynamics at small 
x. This is not surprising. The leading behaviour obtained from BFKL is an x~ x growth, whereas 
for DGLAP we anticipate an increase of the double logarithmic form exp ^A[ln(t/to) ln(l/x)]^ 
where t = ln(Q 2 /A 2 ). However, these are asymptotic predictions. For instance subleading 
In 1/x effects will weaken the BFKL growth in the HERA regime Moreover the DGLAP 
behaviour is dependent on the choice of a non-perturbative input form at some scale Q 2 = Q 2 . 
It has been realized however, that the intimate relation between the increase of the cross 
sections with decreasing x and the absence of transverse momentum ordering, which is the 
basic property of the BFKL dynamics, should reflect itself in the properties of the final states 
in deep-inelastic lepton scattering. Indeed several dedicated measurements have been proposed 
and are being experimentally studied at HERA (see, for example, the reviews in ref. |§). 

The purpose of this paper is to study the detailed properties of the partonic final state 
produced by the gluon emissions along the BFKL chain. In this way we will gain insight into 
the BFKL equation, as well as detailing observables with which to probe the underlying small x 
dynamics. In particular we calculate the decomposition of the (total) deep-inelastic cross section 
into components <J n (n) which correspond to the production of a fixed number n of gluon jets 
each with transverse momentum > fi- That is we study the possible jet configurations in 
the central region between the current jet and the proton remnants. An interesting feature of 
BFKL dynamics is the possibility of producing jets even for \i > Q. One of our aims is to 
quantify the yield of such jet configurations. 
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In the BFKL equation there is a delicate cancellation between the real gluon emissions and 
the virtual contributions. Clearly the cancellation is affected by the resolution > // that 
we impose. In particular we must ensure that the appropriate cancellation between the virtual 
contributions and the "unresolved" real gluon emissions with < fi is maintained throughout 
the calculation. We must therefore first derive a modified form of the BFKL equation which 
will enable us to quantify the number of energetic resolved jets produced along the gluon chain, 
but in which the virtual and unresolved contributions are treated on an equal footing and are 
resummed. This is the subject of section 2. In section 3 we give an analytic solution for the 
resummation at low jet resolution, whereas in section 4 we consider more realistic values of the 
resolution \x and solve the modified BFKL equation by iteration to illustrate the jet decom- 
position of the BFKL gluon. At this stage it is still a theoretical study. In section 5 we use 
the BFKL gluon and the Ay-factorization theorem |KJ to predict the jet decomposition of the 
observable structure function F 2 and the deep-inelastic cross section. Section 6 contains our 
conclusions. 



2. The BFKL equation incorporating jet resolution qx > fi 

In the small x regime the dominant parton is the gluon. Since we no longer have strong- 
ordering in transverse momenta along the gluon chain in Fig. 1 we must work in terms of the 
gluon distribution /(a;, fc T ) unintegrated over its transverse momentum fey. The relation of 
unintegrated distribution / to the conventional gluon distribution is 



g(x,Q 2 



Q 2 dB 



7,2 fi X i ^t) 



(2) 



The unintegrated density / satisfies the BFKL equation which effectively sums up the leading 
as In 1/x contributions. In integral form it may be written [11], [12[ 



= r\v,k 2 T ) + as /" dy> 



d 2 q T 



f(y',k 2 )e(k 2 



4) 



(3) 



with as = Sas/n. We have chosen to use the rapidity variable y = In 1/x instead of x, and so 
the integral in (|j) has come from the replacement 



1 dx' 
x' 



-T - dy'. 



For convenience we have also introduced 



k'^ = \q T + fey | 



(4) 



Note that the dependence on makes the angular integration in d 2 qx non-trivial. The inho- 
mogeneous contribution f(°* in (|3|) corresponds to the "no-rung" contribution of Fig. 1. It is 



2 



the driving term of the equation and has to be input. We implicitly include under the d 2 q T 
integral in (^) the product of theta functions 

0(Q 2 - k'S) Q(k'S - Ql) (5) 

so that the emitted gluon is constrained to the domain Qq < k'S < Q 2 . In the numerical 
predictions shown below we take Qq = 1 GeV 2 and Q 2 = 10 4 GeV 2 . 

Jet structure is embodied in the BFKL equation via real gluon emission from the gluon 
chain prior to its interaction with the photon probe (which takes place through the usual fusion 
subprocess 7$ — > qq). An observed jet is defined by a resolution parameter /i which specifies 
the minimum transverse momentum that must be carried by the emitted gluon for it to be 
detected. For realistic observed jets in the experiments at HERA, the lowest reasonable choice 
for the resolution cut-off parameter /1 appears to be about \x = 3.5 GeV. However, we also 
present results for /1 = 6 GeV and, so as to gain theoretical insight, for the low values of /1 = 1 
and 2 GeV. 

If an emitted gluon has transverse momentum < fx then the radiation is said to be 
unresolved. The unresolved radiation must be treated at the same level as the virtual corrections 
to ensure that the singularities as q\ — > cancel in the q\ integration. To do this we first rewrite 
the BFKL equation (Q) in the symbolic form 

/ = / (0) + dy'K®f{y') } (6) 

J 

where ® denotes the convolution over q?. We divide the real gluon emission contribution into 
resolved and unresolved parts using the identity 

0(4 -/i 2 ) + e(/i 2 -4) = 1, (7) 

where the first term denotes the real resolved emission and the second the real unresolved 
emission. We then combine the unresolved component with the virtual contribution That 
is 

/ = / (0) + f dy'(K R + K uv ) ® f(y'), (8) 
Jo 

where the kernel Kr for the resolved emissions with qx > \x is given by 



K R ® f(y') = a s {k 2 T ) k 2 T I d \ Q(q 2 T - /x 2 ) i f(y>, k'S), (9) 



d 2 q T 

while Kuv, the combined unresolved and virtual part of the kernel, satisfies 



K uv ® f(y') = a s (k 2 T ) J 



d 2 q T 
■nq\ 



2 



k 

T 



-| f(y', k'S) 6(/i 2 - q 2 ) - f(y', k 2 T ) ®{k 2 T - q 2 ) 

rCn 



(10) 



with k'S = \q T + kfl 2 ■ These identifications of the kernels follow by comparing (J8|) with @. 
The q\ — > singularity is now cancelled between the unresolved and virtual contributions, and 
by working with the combined kernel Kuv we will ensure that the cancellation remains intact. 
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We seek a BFKL equation for the real resolved emissions in which the unresolved and 
virtual contributions have been resummed. To do this we write the BFKL equation (j8|) in the 
differential form 

% = (jjf + Kr® f) + K uv ® /, (11) 

and treat the expression in brackets as the inhomogeneous contribution. We solve the inhomo- 
geneous equation in the standard way. We first find a solution to the homogeneous equation 
and then we obtain the full solution via an integrating factor. The homogeneous version of (|TT|) 
is 

dA 

— = K uv ® A (12) 

with solution 

A(y) = exp{yKuv), (13) 
and so the integrating factor is A -1 = exp(— y K uv ). Hence the full solution of ( ]TT| ) is 



f(y) = \l dy'A{y) ® A-\y>) ® + K R ® f(y>) 



V dy'e^ K w ® (y^- + K R ® f{y')^ . (14) 



Thus we have derived a BFKL equation for the gluon distribution / in which the unresolved 
and virtual terms have been resummed in the exponential factor. The equation is of the form 

f(y) = f m (y) + I" dy'K ® f(y') (15) 
Jo 

where the driving term has become 

f {0) (y) = jT dy'e^-y'^ ® ^ (16) 

and the new kernel 

K = e iy - y ' )Kuv ® K R . (17) 

Recall that the original BFKL kernel, K R + Kjjv, has no y (i.e. x) dependence. However, upon 
the resummation of the unresolved and virtual radiation we generate an explicit y dependence. 
In fact the kernel K of ( |T5| ) is a function of only the difference y — y' (i.e. of mar' 'fx) and not y 
and y' individually, see flTTj). 



3. Analytical solution at low fi 

In section 4 we numerically solve the modified BFKL equation for f(y, fcf,) and, by iteration, 
determine the probability of the emission of n gluon jets with qx > fi. However, first it is 
informative to derive an approximate form of the above equation which holds in the (theoretical) 
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limit of small /i 2 / k\. In this limit it is possible to resum the unresolved and virtual contributions 
in a closed analytic form. The crucial observation is that for small /z 2 fk\ we may write 



\q T + k T \ 



in the integrand for the unresolved real emission term in ([H]). Then (fL0|) simplifies to become 



K uv ®f(y') = as(k 2 T )f(y') 



dq\ 
Qt 



e(/i 2 -^) - 9(4 -q 2 ) 



o 



kr, 



-a s {k 2 T ) In ( ^ ] f(y') + O ( | 



(18) 



Thus the homogeneous solution of the BFKL equation (O) is 

A(y) = exp(yK uv ) = exp(-ya s (A£) In^/fi 2 ] 



(19) 



that is the resummation is given by a simple analytic form. As a consequence, in the small /i 
limit, the modified BFKL equation fll~5|) becomes 



f(V,k 2 T ) = fW(y,k 2 T ) + a s {k 2 T ) f V dy'A(y-y',k 2 T ) f ^ e(q 2 - fi 2 ) ^ f(y' , k'S) (20) 



where here k'^ = \ q T + \ 2 , and the driving term is given by 



f {0) (y,k, 



T 



(21) 



Of course for the results presented below we do not use the low /i approximation, although to 
gain insight we will compare the full prediction of ( [TED for /(°) with the approximate 0(fx 2 /k^) 
result given in (|2"T|). 



4. Jet decomposition of the BFKL gluon 

The BFKL equation was expressed in form (|i~5|) specifically so that we can decompose 
the unintegrated gluon distribution / into the sum of contributions with different numbers of 
resolved gluon jets with transverse momenta qr > \i- That is 

oo 

M = E f n (y) (22) 

n=0 

where f n denotes the contribution to the unintegrated gluon distribution / arising from n 
resolved jets in the chain, each with qx > //, see Fig. 2. The n-jet contribution f n obviously 
depends on the resolution whereas the sum / does not. Using fllED we have 

f n (y) = f V dy'K ® P~\y') (23) 
5 



where the O-jet contribution f° = p oS) of (jlB]) and where K is the full resummed kernel of (|I7D . 
For the initial non-perturbative input in ( Jl6| ) we take 

f<®{y) = 3N(1 - e~yf exp(-A; 2 /Q 2 ) (24) 

where the normalization N is fixed so that the gluon, integrated over the region fcf > Qq, 
carries half the momentum of the proton. We set Qq — 1 GeV 2 . 

Although the sum f(y) of (|22|) is independent of /i, the individual contributions f n (y) are 
\x dependent. Recall that ® stands for an integration over d 2 qx (see (|9|) and fllCf)), and that / 
is a function of k\ as well as y. In Figs. 3, 4 and 5 we show the decomposition of f(y, k^) for 
k? = 2, 5 and 10 GeV respectively, in each case taking three different values for the resolution, 
namely ji — 1, 2 and 3.5 GeV. The gluon density, and its decomposition, are not observable 
directly. The choices we have made for /i are, at this stage, solely to gain insight into the 
structure of the BFKL gluon. The results show the following features: 

(i) Gluon jets with /i > k T occur; the probability increases as x decreases. 

(ii) The lower the value of fi, the greater the number of resolved jets, that is the greater the 
preponderance of multijet configurations. 

(iii) As x decreases, the greater the diffusion in lnq^ so that an n-jet configuration first 
increases in probability and then decreases as higher jet-configurations take over. 

(iv) The higher the value of k\ the sooner in x (as x decreases) will a given multijet configu- 
ration go through this rise and fall. 

(v) As fcf.//i 2 increases the 0-jet contribution drops rapidly to zero. 

The results for low values of the resolution parameter [i show that the functions f n have a 
maximum which shifts to smaller values of x with increasing n. This maximum is a straight- 
forward consequence of virtual corrections which, for low fi, are not entirely compensated by 
(unresolved) real radiation. The maximum disappears for large /i and we have this structure 
for all values of fcr- 

Some insight into the behaviour can be obtained from the analytic form presented in section 
3, which applies when /z 2 jk\ is small. In this limit the virtual and unresolved real terms lead 
to a suppression factor 

A(y) = e~ A y (25) 
where A = as ln(k^/fi 2 ). Thus from (|2T| ) we obtain the 0-jet contribution 

f = f(P)f y k l) = e~ Ay f V dy' e Ay ' 3Ne~ k T/ Q o -f- (1 - e~ y 'f , (26) 

Jo ay' 

that is the kj- dependence of f° is essentially the same as the kj- dependence of the driving 
term f(°> of (p4|). This explains the origin of feature (v), that the 0-jet contribution falls 
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rapidly to zero with increasing k\. 



Fig. 6 compares the analytic approximation with the full 
result for /z 2 = 1 GeV 2 and k\ = 4 GeV 2 . We see that the analytic form reproduces the 
shape of the numerical solution, but fails in the normalization. Also the peak in the numerical 
prediction shifts slightly to smaller x. Thus the analytical approximation cannot be used as a 
valid representation for the jet contributions, even for a resolution as low as 1 GeV 2 . 



5. Jet decomposition of Fi(x,Q 2 ) at small x 

We are now in a position to estimate the probability of the different multijet configurations 
in the small x observables that are driven by the BFKL gluon. The most relevant process to 
study is deep-inelastic scattering at HERA. Using the results of section 4, we calculate the jet 
decomposition of the proton structure functions Fi(x,Q 2 ). In other words we determine what 
fraction of events that make up the inclusive measurement of Fi(x, Q 2 ) contain no-jets, one jet, 
two jets etc. as a function of x, Q 2 and the jet resolution parameter fi. Recall that our jets are 
gluons emitted with transverse momentum > /i. 

From knowledge of the BFKL gluon / we can determine the behaviour of the structure 
functions via the kr factorization theorem, see Fig. 2. For the transverse and longitudinal 
functions we have 



T,L 



x,Q 2 



1 dx' 



x 



f 



x 



X 



I L.2 ^ 



Q 1 



(27) 



where, to lowest order, photon-gluon fusion F 19 is given by the quark box (and crossed box) 
contributions, as shown in Fig. 7. To carry out the integration over the quark line in Fig. 7 we 
express its four momenta k in terms of the Sudakov variables 



ap 



(3q + n T 



where q' = q + xp and p are the basic light-like momenta (q and p are the 4-momenta of 
the virtual photon and proton respectively). The variable a is fixed by the quark mass-shall 
constraint, leaving integrations over (3 and Kt- Evaluating the box contributions, equation fl27D 



then becomes 14 



F T (x,Q 2 



Fl(x,Q 2 ) 



2E 



Q 1 



„ q 47T Jkl If 



J q d/3 J d 2 n' T a s f Qj,A;, 




m„ 



D X D 7 



(28) 



(29) 
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where the denominators 



D 1 = 4 + 0(1 - 0)Q 2 + m 2 

D 2 = |kt- k T \ 2 + /5(1 - /?)Q 2 + m. 



and where /t^ = — (1 — (3)k,T- The x' integration of (^71) is implicit in the d 2 n' T and 
integrations. Indeed x' is fixed in terms of n' T and 



/ 

x 



k'j. + m 2 k T 1 ' 



1 H f h — 

+ - 0)q 2 + q 2 



(30) 

which ensures that the requirement < x' < 1 is satisfied. Of course the integration regions of 



(28) and (p9| ) must be constrained by the condition 

x'(P,K%,k T ,Q 2 ) > x (31) 

so that the argument z — x/x' of / satisfies the requirement z < 1. In ( f28|) and ( p9|) we sum 
over the quark flavours; we take the masses to be m q = for u, d, s quarks and m c — 1.5 GeV 
for the charm quark. The argument of as is taken to be + Wq, which allows integration 
over the entire region of k'^. For the light quarks we take mo = 1 GeV 2 ; the results are not 
very sensitive to variations of mo about this value. For the charm quark contribution we set 



mr, = mz. Also we set = 1 GeV 2 . 



The jet decomposition of Fl,t are simply obtained by substituting the n-jet unintegrated 
distribution f n into ( p8l) and (p9|) . In this way we can break down the observables into their 
component n-jet contributions, for example for F 2 = + Ft we have 

oo 

F2 = E F 2- (32) 

n=0 

Figs. 8 and 9 show the components F£(x, Q 2 ) for deep inelastic events containing n observed 
jets, where in the upper plot we require the jets to have qr > 3.5 GeV, whereas in the lower 
plot we demand qr > 6 GeV. Figs. 8 and 9 correspond to Q 2 = 10 and 20 GeV 2 respectively. 
For these choices of jet resolution it can be seen that the 0-jet configuration dominates. That 
is most of the emission from the BFKL ladder is in the form of unresolved and virtual gluon 
radiation. As expected the n-jet configurations first become important (with decreasing x) for 
the lower resolution, /i = 3.5 GeV, and the higher Q 2 value, Q 2 = 20 GeV 2 , and begin to 
compete with the 0-jet rate for x <• 10~ 5 . In fact the 4-jet rate becomes comparable with the 
0-jet rate for x ~ 10 -6 . 

Although the 0-jet configuration dominates in the HERA kinematic regime, there is still a 
non-negligible contribution from resolved jets. For example, at Q 2 = 10 GeV 2 and x = 2 x 10~ 4 
the 1 and 2 jet contributions are each approximately | of the 0-jet rate, and even the 3 and 4 
jet configurations occur at a reasonable rate. Also notice the production of resolvable jets with 
fJ 2 ~ Q 2 is important — this is a straightforward consequence of diffusion in k T . 
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The experiments at HERA show that the (inclusive) structure function F 2 rises as x de- 
creases. How is this rise made up from the various multijet configurations? First we look at 
the results for the lower jet resolution, /i = 3.5 GeV. Although the 0-jet rate dominates, its 
increase with decreasing x is relatively weak compared to the data. The rise of F 2 comes from 
the increasing importance of the higher jet configurations. On the other hand at the higher 
resolution, /i = 6 GeV, the 0-jet configuration is even more dominant and shows a steeper rise 
over the same x range, as is required for consistency of the results. This characteristic difference 
should hopefully be seen in the measurement of the individual jet structure functions. 

The cross section for deep-inelastic scattering is readily calculated from F TL . We have 

o = W J ^ J ^ {y'xF^Q 2 ) + (l-y)F 2 (x,Q 2 )} (33) 

where as usual y = Q 2 /xs, F T = 2xFi and F L = F 2 — 2xFi. Here we present results for the 
component cross sections a n for deep-inelastic events containing n jets with q T > /i, again for 
two choices of resolution /i = 3.5 and 6 GeV. We take = 300 GeV and integrate a over the 
interval 0.01 < y < 0.5 so as to approximately reproduce the HERA domain. Figs. 10, 11, 12 
and 13 show respectively the 0, 1, 2 and 3 jet cross sections integrated over x and Q 2 bins of 
size Ax = 2 x 10 -4 and AQ 2 = 10 GeV 2 , where the two entries in each bin correspond to a 
gluon jet with resolution fj, = 3.5 and 6 GeV respectively. We see that there are an appreciable 
number of identifiable jets. For example, if we take a resolved jet to be one with > 3.5 GeV 
and an integrated luminosity C = 10 pb^ 1 , then in the bin defined by 0.8 x 10~ 3 < x < 10~ 3 
and 15 < Q 2 < 25 GeV 2 we predict 1052, 790, 392 events containing 1, 2, 3 jets as compared 
to 6790 events with no identifiable jet. 

Recall that the predictions are obtained by numerically solving the BFKL equation for the 
gluon. The normalisation is dependent of the choice of the cut-off. Here we have taken the 
cut-off to be 1 GeV 2 , which was found to give a satisfactory description of the inclusive F 2 
distribution. However, the fraction of events containing 0, 1, 2, . . . identifiable gluon jets is 
independent of the choice of the cut-off. For example, for the above (Ax, AQ 2 ) bin and for the 
lower jet resolution of fi = 3.5 GeV we find 75% of the cross section contains no observable jet 
and that 1, 2 and 3 jets occur 11, 8, 4% of the time respectively. Only 2% of the events contain 
more than 3 jets. For the higher jet resolution of /i = 6 GeV we predict that the BFKL chain 
will give 90% of the events with no observable jet, leaving only 10% of the total to be split 
between 1, 2, . . .jet events. 

We see from Figs. 11 and 12 that the 2-jet rate is comparable to the 1-jet rate, and more- 
over that the 2-jet/l-jet ratio increases with increasing resolution /i. This type of behaviour is 
consistent with the expectations of the conservation of transverse momentum. 

6. Summary and Conclusions 
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In this paper we have formulated a modified form of the BFKL equation which allows an 
exclusive analysis of the multijet yields in deep-inelastic lepton scattering in the small x regime. 
The jets are defined as gluon emissions from the BFKL chain which have transverse momenta qx 
greater than a specified resolution /i. We first solved the modified BFKL equation to determine 
the jet decomposition of the unintegrated gluon distribution f(x,k^). We then used the fc*r- 
factorization theorem to determine the jet decomposition of the structure function F 2 (x, Q 2 ) 
and of the total deep-inelastic cross section in the HERA small x regime. We presented the 
jet decompositions as a function of the kinematic variables and for different choices of the jet 
resolution parameter \i. 

The modified BFKL equation is shown symbolically in ( |i5| ) and the kernel K in ([17]). The 
equation embodies a resummation of the virtual contributions together with the unresolved 
real gluon emissions with < fi. As a consequence the kernel K has an explicit y = lnl/x 
dependence, which depends on the amount of unresolved radiation and so is a function of \x. 
Indeed for unrealistically low values of /i we derived, for pedagogic purposes, the analytic form 



of the y dependence of the kernel, see (19). For the more realistic numerical solutions that we 
present the correlation between the x dependence of the n-jet cross sections and the resolution 
parameter /i is apparent. 

The behaviour of the n-jet contribution to the gluon /, or to F 2 , exhibits a characteristic 
behaviour as x decreases, rising to a maximum and then falling back to zero. The higher the 
value of n the lower the value of x at which the maximum occurs. In the HERA small x 
regime the behaviour is only apparent for low choices of the parameter fi, for example /i ~ 1 
GeV, see Figs. 3, 4 and 5. For experimentally realistic values of the resolution parameter (say 
fi = 3.5 or 6 GeV) the maxima shift to very small values of x. The dominant contribution in 
the HERA range then comes from events with no resolved gluon jets emitted from the BFKL 
chain. Nevertheless the 1, 2, 3, . . .jet rates are still significant. An interesting feature of the 
multijet cross sections is that they are non- negligible even if fi > Q. The existence of such jets 
with qx > Q is a straightforward consequence of the characteristic Ink 2 -, diffusion along the 
BFKL gluon chain. 

To sum up, we have made an exploratory study of a form of the BFKL equation which 
allows the final state jet configurations to be determined in a consistent manner. We solved 
the equation and presented sample results to illustrate the properties of these gluon jets which 
occur in deep inelastic scattering at small x as a result of BFKL dynamics. Of course there are 
subleading lnl/x and fixed-order QCD jet contributions to consider. These may modify the 
predictions in the HERA regime, but with decreasing x the BFKL behaviour should become 
increasingly dominant. One non-leading effect is the imposition of the constraint q 2 T < x n k\jx 
(in the notation of Fig. 1) which follows from the requirement that the virtuality of the gluon 
links is dominated by —k\ 0, §. If this were done we find that it would limit the available 
phase space for multijet production and, as a consequence, reduce the yield of multijet events. 
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Appendix: Numerical techniques used to solve the BFKL equation 

Here we briefly describe the numerical method that we used to solve the BFKL equation 
(15) . The starting point is the Chebyshev polynomial expansion of the unintegrated gluon 
distribution f(y,k^) in which we map the region Ql<k^< Q"j into the interval (—1,1) in 
terms of the variable r defined by 



r(4) = 21n('JU /J% 



\ QfQo J I \Q 

We expand the gluon distribution / in the polynomial form 



(Al) 



N 

f(v,k 2 T ) = £ a(r(4)) My) (A2) 



1=1 

where fi(y) are the values of f(y, k%>) at the {kj)i nodes obtained from 

(gjj = (QA n 
QfQo \Qo) ' 

with Tj defined by 



(A3) 



Ti = cos [(, _ l)n/N], (A4) 

and the number of terms in the Chebyshev polynomial. The k\ dependent functions Cj are 
obtained from the Chebyshev polynomial functions 



T n {r) = cos|^n arccos(r)J (A5) 

and are given by 

2 N 

Ci{T) = E VnTn^T^), (A6) 

iV n=l 

where u n = 1 for n > 1, and V\ = A good approximation for the k\ dependence of / is 
obtained with typically A^ = 20. 

The expansion ( |A2| ) is then substituted into the BFKL equation ( |T5D to give the discretised 
(symbolic) form 

rv N 

My) = f 0) (y) + dy'Y, K hk (y-y')fk(y'), (A?) 

11 



where the full kernel (jlTf) now becomes 

Ki,k = E [e {y - y ' )Kuv h K? k (A8) 
i 

and the input distribution fl (y) of (]X6|) is 

/^(v) = [ dy> E [e^'^k* (A9) 

The substitution of (|A2|) into (|7|) and (|[) gives the explicit form of the kernels K R and -fQ/v 
respectively. The BFKL equation (|A7|) is a Volterra-type integral equation, which we solve 
iteratively for the fi(y)'s. The gluon distribution f(y, fcf,) is then reconstructed from ( |A2|) . 

We also use a Chebyshev interpolation to calculate the Y = y — y' dependence of the 
matrix elements of exponential matrix in (|A8|) and (A~9). For convenience we denote the matrix 
elements 

[e YKuv ] i>k = M(Y) i>k . (A10) 
As before we expand in terms of Chebyshev polynomials 

M{Y) hk = E CMY)) Hk ( An ) 

where M J are the values of M(Y) at the nodes Yj. Here we take J = 10. We map the relevant 
region < Y < F max , where Y max = log(l/X min ), into the interval — 1 < r < 1 by choosing 

t(Y) = (2Y - y max )/y max . (A12) 



The Cj is given by ( |A6| ) (together with (|A4|) and (|A5|)) with i replaced by j. It remains to 
calculate M(Y) at the nodes Y — Yj. We do this by solving 

dJ ^ 1 = >>rr)„ ; U,MY). (A13) 
using the Runge-Kutta method with the boundary condition M ijk (Y = 0) = i^. 
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Figure Captions 

Fig. 1 The unintegrated gluon distribution, f(x, k^), is effectively the sum of the ladder diagrams 
formed by the modulus squared of such amplitudes. The leading a^lnl/x resummation 
is accomplished by the BFKL equation. 

Fig. 2 The modulus squared of this diagram gives the component F™ of the proton structure 
function Fi which arises from the contribution f n to the gluon distribution / in which 
there are n resolved gluon jets emitted along the BFKL chain, that is n gluons with 
qr > fi. The black circles are to indicate the presence of both virtual and unresolved 
gluon emissions. The component F™ is calculated by the /^-factorization theorem, which 
has the symbolic form F t n = F? 9 <g> f n , see (E7I) and (R2p . 

Fig. 3 The n-jet contributions to the unintegrated gluon distribution f(x, kj) for three different 
values of the jet resolution parameter fi and for kr = 2 GeV. 

Fig. 4 The same as Fig. 2 but for fcy = 5 GeV. 

Fig. 5 The same as Fig. 2 but for = 10 GeV. 

Fig. 6 The comparison of the analytic and numerical solutions for the 0-jet contribution f°(x, kj) 
to the unintegrated gluon distribution. 

Fig. 7 The quark box and crossed box diagrams describing photon-gluon fusion, F 19 in (P7|). 

Fig. 8 The decomposition of the proton structure function F2(x,Q 2 ) into contributions coming 
from different numbers of resolved gluon jets for experimentally accessible values of the 
resolution parameter fi = 3.5 and 6 GeV. The decomposition is shown as a function of x 
for Q 2 = 10 GeV 2 . 

Fig. 9 The same as Fig. 8 but for Q 2 = 20 GeV 2 . 

Fig. 10 The cross-section (in pb) for deep-inelastic scattering in which there are no resolved gluon 
jets shown in different x, Q 2 bins in the region accessible at HERA. The width of the bins 
are AQ 2 = 10 GeV 2 and Ax = 2 x 10~ 4 . The upper and lower values correspond to the 
resolution parameter fi = 3.5 and 6 GeV respectively. 

Fig. 11 The same as Fig. 10 but from the contribution in which there is one, and only one, gluon 
resolved jet with q T > fi. 

Fig. 12 The same as Fig. 10 but for two resolved gluon jets. 

Fig. 13 The same as Fig. 10 but for three resolved gluon jets. 
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